Effective Temperatures in Driven Systems: Static vs. Time-Dependent Relations 
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Using simulations of glassy systems under steady-state shear, we compare effective temperatures 
obtained from static linear response with those from time-dependent fluctuation-dissipation rela- 
tions. Although these two definitions are not expected to agree, we show that they yield the same 
answer over two and a half decades of effective temperature. This suggests that a more complete 
conceptual framework is necessary for effective temperatures in steady-state driven systems. 
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Temperature is one of the fundamental variables in an 
equilibrium system that determines not only the system's 
average properties, such as pressure or density, but also 
fluctuations around those averages. Temperature also 
relates, via linear response, fluctuations in a thermody- 
namic quantity to that quantity's response to a small 
perturbation in its conjugate variable. When a system 
is far out of equilibrium, temperature is no longer well- 
defined. Nevertheless, in many cases there will still be 
fluctuations although they are not thermal in origin. An 
example of this is a steady-state driven system such as 
a sheared material where shear introduces fluctuations 
that are not described by a thermal bath temperature. 
Can one define an appropriate "effective temperature" to 
characterize these fluctuations? 

For the idea of an effective temperature to be use- 
ful, a clear prescription for defining it should exist and 
this prescription should apply generally in different con- 
texts. Various groups have defined and measured dif- 
ferent effective temperatures in systems far from equilib- 
rium Ql QUI Q| . We will show here that two prescriptions 
based on linear response that have seemed to be fun- 
damentally incompatible do, surprisingly, give the same 
value for the effective temperature. For this to occur, nei- 
ther prescription can work infallibly in all cases. Thus, 
our results pose the conceptual question: under which 
conditions should either linear response prescription be 
applied? 

Linear response provides a large number of possible 
definitions for the effective temperature, each based on 
a different pair of conjugate variables, which all reduce 
to the true temperature in thermal equilibrium. Effec- 
tive temperatures based on these different relations have 
been used to model simulations of particulate systems 
driven out of equilibrium by steady-state shear [Bl EJ • 
There are two ways in which such relations have been 
implemented with significant success. One way uses 
static linear response which relates equal-time fluctua- 
tions to the response at infinite-times to yield an effec- 



tive temperature Tj. The other way, argued to be more 
fundamental^,!^], measures the autocorrelation function 
and relates it to the response as a function of time. In 
equilibrium there is a strict proportionality between cor- 
relation and response at all times and therefore a single 
well-defined temperature, but in driven systems the situ- 
ation is more subtle. The conceptual picture behind this 
is that there can be two widely separated time scales in 
the presence of shear. Degrees of freedom that decay on 
a short time scale are characterized by one effective tem- 
perature, Ts, while those that take a longer time to decay 
are characterized by a different (higher) value, Tj,. It has 
been predicted that Ts is the bath temperature because 
fast degrees of freedom decay before shear has any effect, 
while Tl , which can only be obtained from the long-time 
behavior of correlation and response, corresponds to a 
well-defined effective temperature characterizing struc- 
tural rearrangements driven by shear 

As will be made clear below, an effective tempera- 
ture based on infinite time, or static linear response, 
Tf, should generally not agree with Tl obtained from 
time-dependent linear response in non-equilibrium sys- 
tems. Indeed, for the same pair of conjugate variables 
we always find in our simulations that Tj =/= Tl- How- 
ever, for a wide variety of simulations we find that Tj 
for one conjugate pair can be equal to Tl for a different 
pair. Therefore, it appears that under different condi- 
tions both definitions of effective temperature must be 
equally valid. The conceptual framework that has 

been used to argue for the validity of Tl would imply 
that Ti should never be valid, in contradiction to our 
findings. That scenario must therefore be incomplete. 

In order to demonstrate these results, we have per- 
formed numerical simulations of systems undergoing lin- 
ear shear flow in both two and three spatial dimensions 
(2d and 3d) under conditions of fixed volume, fixed num- 
ber of particles and fixed shear rate. The systems are 
composed of 50-50 bidisperse mixtures with diameter ra- 
tio 1.4, which prevents crystallization and segregation. 
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FIG. 1: (a) R p ^)(t) vs - correlation C p n,)(t) in a 2d system 
with repulsive harmonic interactions at a bath temperature 
Ta_b = 10~ 4 , shear rate 7 = 0.01, and packing fraction cj> — 
0.90 for N = 256 particles. The wavevector k lies in the 
shear gradient direction; k — 4.5, 9, and 11 are shown with k 
decreasing from top to bottom. The solid lines are guides to 
the eye and have slopes equal to — 1/Tl- (b) Rp(t) vs. Cp(t) 
for pressure in the same system as (a). The solid line is a 
guide to the eye. It has slope — 1/Tl and intercept 1/Tl 
where Tl is obtained from (a) . Note that Ti from (b) is equal 
to Tl from (a). The results are in the linear response regime, 
as shown by the open circles and dashed lines in each of the 
figures, which are for two magnitudes of the perturbing field 
that differ by a factor of 5. 

The system is enclosed in a cubic simulation cell with 
Lees-Edwards periodic boundary conditions to impose 
shear in the ir-direction and a shear gradient in the y- 
direction. Particles interact via one of the following pair- 
wise, finite-range, purely repulsive potentials: 



V hs {r l3 ) 



V RLJ {n 3 ) 



-(l-rij/aij) 2 
2c 



(l-nj/aij) 



5/2 



(1) 



72 [K-A y ) 12 -2(a y /r 2 ,) 6 + l] 



where e is the characteristic energy scale of the interac- 
tion, Uij is the average diameter of particles i and j, and 
Tij is their separation. All potentials (harmonic spring, 



respectively where m is the particle mass and a is the 
small-particle diameter. 

We have studied both thermal and athermal (or dis- 
sipative) systems to show that our results are not spe- 
cific to any particular dynamics. Thermal systems under 
shear can be described by the Sllod equations of motion 
for the position fi and velocity fluctuation u, of each par- 
ticle around the average linear velocity profile 12]: 



~dt 



~dt 
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where F[ — ~J2j^( r ij)/^ r ij^ij ^ s the repulsive force 
on particle i due to neighboring particles j, 7 is the shear 
rate, and a is chosen to fix the kinetic energy per degree 
of freedom, Tke, associated with velocity fluctuations. 
We always set Tk e to be below the glass transition tem- 
perature of the unsheared system. 

Athermal dissipative systems can be described bv[l3|: 
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where v\ is the total velocity (including shear) of parti- 
cle i, b > is the damping coefficient, and the sum over 
j only includes particles that overlap particle i. At fi- 
nite shear rate, these systems reach a steady-state where 
the power put in by the shear flow balances the power 
dissipated. In this study, we focused on underdampcd 
dissipative dynamics and therefore fixed the dimension- 
less damping coefficient b* = ba/^em ^ l. 

We now describe how to calculate the effective tem- 
perature from a set of conjugate variables and then show 
explicitly the incompatibility of the static and the time- 
dependent definitions for the same conjugate pair of vari- 
ables. Since we are concerned with systems in steady- 
state shear, we assume a final steady state in which av- 
eraged quantities become time independent. Consider 
an observable, A(t), that fluctuates in time t, such as the 
number density or the total pressure of the system. Then 
one can define the autocorrelation function: 



C A (t) = (A(t)A(0)} - (A)' 



(4) 



where (A) represents an average of A over time and con- 
figurations. If B is the thermodynamic field conjugate to 
A, then one can also define an integrated response func- 
tion which measures the response to a small constant 
perturbation, SB, applied from time t = onward: 



R A {t) 



(A(t)) - (A(0)) 
SB 



Hertzian nonlinear spring and repulsive Lennard- Jones) If wg nQW introduce the rescaled variab i es 



are zero when ry > c.y . Our results have been obtained 
with packing fractions ranging from cf) = [0.70, 0.84] in 
3d and cf> — [0.85, 1.20] in 2d, which are all above ran- 
dom close-packing ^l| . We varied the number of parti- 
cles in the range N = [256, 1024] and found no apprecia- 
ble finite-size effects for the results reported here. The 
units of length, energy, and time are a, e, and <j\J m/e, 
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the fluctuation-dissipation relation states that 
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FIG. 2: (a) i? p ( fc )(t) plotted versus correlation C p ^){t) in a 
2d sheared, athermal system with repulsive harmonic interac- 
tions at cj) = 0.90 and 7 = 0.01. Again, k lies in the shear 
gradient direction and k = 4.5 and 9 are shown with k de- 
creasing from top to bottom. The solid lines have slope equal 
to — 1/Tl. (b) Rp(t) plotted versus the autocorrelation func- 
tion Cp(t) for pressure P in the same system as (a). The 
straight line has slope —1/Tl and intercept 1/Tl where Tl is 
obtained from (a). 

This implies that if Ra is plotted parametrically against 
Cai the result should be a straight line with slope —1/T 
for an equilibrium system. Moreover, the infinite-time 
limit, which is the intercept of such a plot on the R- 
axis, has the value 1/T. According to Eqs. Hj]and[7| the 
intercept satisfies 



R A (t = 00) = C A (t = 0)/T 



(8) 



Thus, the intercept, which defines a temperature Tj, cor- 
responds to the static linear response relation whereby 
the infinite-time response, RA{t — 00) (Eq. |SJ, is related 
to the equal-time correlation function, CU(i = 0) (Eq.0}. 

As indicated earlier, Kurchan |8( has predicted that for 
driven systems in steady state such a parametric plot has 
two regimes 0, Q . Berthier and Barrat have conducted 
simulations of sheared Lennard Jones glasses and shown 
that at short times (Ca close to one and Ra close to zero), 
the slope of the line, —1/Tg, defines a temperature char- 
acterizing the fast modes in the system and corresponds 
to the bath temperature, Tke- At long times the slope, 
— 1/Tl, is a good measure of the effective temperature 
produced by the shear for the slow modes. We show a 
similar plot for a 2d sheared thermal system in Fig. Of a) 
using as the variable the Fourier component of the num- 
ber density of the large particles at various values of the 
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For this variable, we calculate the incoherent part of the 
scattering function for the large particles, C p nA(t) In 
Fig. nja) , there is a well-defined slope at short times and 
a smaller slope at long times (i.e., at smaller values of 
CU(i)). It is quite obvious that the long-time slopes are 
all the same so that there is a common effective temper- 
ature that describes the fluctuation and response for all 
of these variables. These results are completely consis- 
tent with the results found by Berthier and Barrat [f| 
for a three-dimensional system. Note that in all of these 
cases, the value at which each of these curves intercepts 
the Rp(k) axis cannot have the value 1/Tl. This could 
only be the case if Ts = Tl (as in equilibrium), or if the 
regime corresponding to Ts shrinks to zero If the 

curve is not a straight line, then Tj (determined from the 
intercept) must be different from either Ts or Tl- 

In Fig. ^b) , we show the parametric plot for the iden- 
tical system as in Fig. Hfa) but for a different variable, 
namely the total pressure, P. We calculate the pressure 
P = Paa I d in d spatial dimensions, using the following 
expression for the pressure tensor 
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where a,/3 = x,y, or z and L is the edge length of the 
simulation box. The shape of this curve is very different 
from those shown in Fig.QIa). In this case, the response 
rises rapidly at short times and then turns over and be- 
comes horizontal at long times The striking result 
is that, although this curve is manifestly different from 
those shown earlier, it has an intercept temperature, Tj, 
with the same value, Tj = Tl, obtained from the late- 
time slopes of the curves in Fig.^a). 

We stress here that this is not a coincidence. In Fig. [21 
we show for an athermal system the response versus cor- 
relation plots for the same sets of variables as shown for 
thermal systems in Fig.^ Again we see that Tj obtained 
from P has the same value as Tl obtained from p(k). To 
indicate the full extent of agreement we plot in Fig. 
the ratio Tl/Tj versus log 10 Tl for all systems studied. 
This figure shows that within error Tl = Tj over two and 
half decades of effective temperature. These data are col- 
lected from thermal and athermal simulations at different 
values of the shear rate and density in 2d and 3d for sys- 
tems with different particle interactions. For the thermal 
simulations, we also varied the bath temperature. 

We have shown that both static linear response 
and time-dependent fluctuation-dissipation relations can 
yield consistent values of effective temperature. How- 
ever, it is unclear when one should use static relations 
and when one should use time-dependent ones. Pressure 
is not the only observable for which the static relation is 
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FIG. 3: The ratio T h jTi is plotted vs. T L , where Tj is ob- 
tained from the 7?-intercept of integrated response (R) vs. 
correlation (C) for pressure and Tl is obtained from the long- 
time slope of R vs. C for p(k). The largest error comes from 
estimating the long-time slope of R vs. C for p(k). 

appropriate. Previously, we showed that static relations 
yield a consistent effective temperature for shear stress 
and potential energy, as well [7j. It is also not true that 
one should invariably use the static relation for quantities 
involving pressure. We have calculated response vs. cor- 
relation for pressure at different values of k in 2d sheared 
systems. At large k, Tl yields the correct effective tem- 
perature whereas at k — Tj does. 

It is not true that the effective temperature for all 
zero-wavevector observables should be given by static lin- 
ear response. While static linear response appears ap- 
propriate for k = pressure, energy and shear stress, 
we find other examples for which this does not hold, 
including P zz |5j and the deviatoric pressure Pdev — 
(2P U - P 22 - P 33 )/3 where 1,2,3 = x,y,z M- We 
find that the response-correlation curves for each of these 
k = variables are not flat as in Fig. IHb), but have 
nonzero long-time slopes. However, the corresponding 



effective temperatures can be off by factors of 5 — 10 
from those shown in Fig. ^a) even when other effective 
temperatures agree. 

There are regimes where the idea of an effective tem- 
perature is valid and regimes where it breaks down |l7j . 
For example, at high densities (e.g. <j> w 1.1, typical of 
liquids) and high bath temperatures (T > 0.1T g ) we find 
that all the different effective temperatures are the same. 
However, these begin to deviate as the bath tempera- 
ture is lowered. Previous studies H, Q have suggested 
that the concept of effective temperature should be valid 
when there is a clear separation between the short-time 
(bath temperature) regime and the long-time, shear-rate- 
dependent regime. This criterion cannot be sufficient be- 
cause as the bath temperature is lowered, the separation 
between these two time scales does not decrease. 

We have shown that Tj for the zero-wavevector pres- 
sure is equal to Tl for p{k) over a range of two and a half 
decades in effective temperature. This remarkable result 
suggests that static as well as time-dependent linear re- 
sponse relations can be used to define a consistent effec- 
tive temperature, in contradiction to expectations based 
on spin models subjected to non-conserved fields @, El- 
These results also leave us with a puzzle: when should 
one use static linear response and when should one use 
a time-dependent relation? For a given pair of conjugate 
variables, there is no obvious criterion for which of these 
two kinds of relations should be used. 
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